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Abstract 

Edge-to-edge interprecipitate distance distributions are critical for predicting precipitation 
strengthening of alloys and other physical phenomena. A method to calculate this 3D distance 
and the 2D interplanar distance from atom-probe tomographic data is presented. It is applied 
to nanometer-sized Cu-rich precipitates in an Fe-1.7 at.% Cu alloy. Experimental interprecipitate 
distance distributions are discussed. 
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Many physical properties of materials depend on the edge-to-edge interprecipitate dis- 
tance, A e _ e . The applied stress required for a dislocation to glide past or climb over pre- 
cipitates depends on A e _ e [1], as does precipitate coarsening and electrical conductivity. 
Frequently, A e _ e is merely approximated by assuming the precipitates form a cubic array or 
a square array in a plane [2]. It is also assumed that precipitates are spherical with a known 
precipitate size distribution (PSD) (usually either all precipitates are the same size or they 
obey the PSD derived by Lifshitz and Slyozov [3] and Wagner [1] (LSW)). Real materials 
are almost always more complicated. 

Much of the past work on calculating the distance between precipitates or other mi- 
crostructural features of interest [5, 6] (whether interprecipitate distances [2], mean free 
paths or chord lengths [7], or nearest-neighbor distribution functions [8, 9, 10, 11]) has 
been theoretical. Experimental characterization of A e _ e requires a microscopic technique 
that has: (i) a high enough spatial resolution to define clearly each and every precipitate; 
(ii) a large enough analysis volume to capture many precipitates and to exclude boundary 
effects; and (iii) 3-dimensional information (without suffering from precipitate overlap or 
truncation). For nanometer-sized precipitates, the local-electrode atom-probe (LEAP®) to- 
mograph (Imago Scientific Instruments) satisfies these requirements [12, 13]. Despite these 
capabilities, it has not been previously utilized to gather this information and the little 
available experimental data for A e _ e comes from 2D techniques. These cannot be com- 
pared directly to models of 3D microstructure, but only to 2D slices from theoretical 3D 
microstructures [14]. 

In this article, an algorithm to calculate A e _ e from LEAP tomographic reconstructions 
is presented and applied to a binary Fe-Cu alloy. This alloy and many other steels are 
strengthened by a high number density of nanometer-sized copper-rich precipitates [15]. 
Many of the proposed precipitate strengthening mechanisms depend on A e _ e [16, 17]. Atom- 
probe tomography has been used to study the size, morphology, and chemical composition 
of Cu precipitates [18, 19, 20, 21, 22], but not to measure A e _ e . An Fe-1.7 at.% Cu alloy was 
solutionized at 1000°C for 1 h and 845°C for 6 h. It was subsequently aged for 2 h at 500°C. 
This treatment leads to a high number density ((1.2 ± 0.1) x 10 24 m -3 ) of nanometer-sized 
precipitates (with a mean radius, (R), equal to 1.3 nm). The specimens were cut, ground, 
and then electropolished into tips. The LEAP tomographic experiment was conducted with 
a 50 K specimen temperature, a 5-10 kV specimen voltage, pulse fraction of 15%, and a 
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FIG. 1: (a) A LEAP® tomographic reconstruction of an Fe-1.7 at.% Cu specimen, whose thermal 
history is detailed in the text. Only Cu atoms are displayed for clarity, (b) The 546 precipitates 
are fitted as ellipsoids [23] (c) A Delaunay mesh connects the precipitate centers. This is used to 
find "interacting precipitates" and to exclude the convex hull. 

pulse repetition rate of 200 kHz to collect ca. 1.3 x 10 6 ions in a 148 x 66 x 62 nm 3 volume 
(Fig. la). The computer program ivas (Imago Scientific Instruments) was used to analyze 
the data. Precipitates are isolated using a modified envelope algorithm [24]. Because Cu 
partitions strongly to precipitates [21], an isoconcentration surface was not necessary to 
distinguish the 546 precipitates in this dataset. 

The interprecipitate distance algorithm begins by representing these precipitates with 
simpler geometric shapes. While A e _ e between spheres is simple (it being the difference of 
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the center-to-center distance and the precipitate radii), spheres do not adequately represent 
many precipitate morphologies. Instead, best-fit ellipsoids to the precipitates are calculated 
(Fig. lb) employing a recently presented algorithm [23]. The 4x4 transformation matrix 
calculated with that algorithm translates, rotates, and scales a unit sphere centered at the 
origin to an ellipsoid that preserves the centroid, principle axes, and moments of inertia of 
a precipitate. 

A Delaunay tedrahedral mesh is generated from the precipitate centroids (Fig. lc) [25, 26]. 
The Delaunay mesh is the geometric dual of the Voronoi diagram; mesh segments connect 
neighboring precipitates whose Voronoi cells touch. It decreases the number of precipitate 
pairs for which A e _ e is calculated to a group of neighbors. The mesh also finds the 75 
precipitates that make up the convex hull. These outer-most precipitates are allowed to 
be nearest neighbors of the inner precipitates, but their own nearest neighbors are not 
calculated, as they might fall outside the volume of the analysis. 

The distance between two ellipsoids is found utilizing the constrained optimization by 
linear approximation (COBYLA) algorithm [27]. This general optimization algorithm is 
chosen over more efficient algorithms that calculate explicitly the distance between ellip- 
soids [28, 29], so that it can be used with other abstractions of precipitate morphology (such 
as the isoconcentration surface) and additional constraints (such as calculation of interplanar 
edge-to-edge distances) and because a gratis implementation exists [30]. COBYLA mini- 
mizes the distance between two points, x and y in the analysis space, y^2j {xj ~ Vjf ■ The 
constraints are that x and y must fall on the ellipsoid. This is simplified by the fact that 
applying the inverse transformation of ellipsoids transforms them back into unit spheres, 
centered at the origin (so • xj x = 1 and . yj v = 1, where the superscript Tj is the 
inverse transform of the best-fit ellipsoid for precipitate i). The initial guess is chosen as 
the two closest points that satisfy these constraints that lie on the line that connects the 
precipitate centers. Interprecipitate distance distributions (IDDs) may be generated using 
different combinations of Delaunay neighbors, as in Fig. 2. An IDD is the convolution of a 
PSD and the center-to-center distances. In Fig. 2a, an IDD for all 6,671 Delaunay neighbor 
distances yields a mean 3D interprecipitate distance, (Ag? e ), of 16 nm. Figure 2b displays 
two subsets of this IDD, each with 471 lengths. The distance between nearest precipitates is 
often used to calculate precipitate-dislocation interactions. The IDD for this is much sharper 
and (Ag^ e )=2.6 nm. Precipitates that are very close to one another might be bypassed as 



4 



0.08 



;>>0.06 



£ 0.04 

o 

0.02 



0.00 
0.25 



0.20 

o.i5 

*8 

o o.io 



0.05 



0.00 



20 



40 



(a) 



J5P 



PI All Dclaunay distances 
(A^ e )=16 run 



(b) 



Shortest Dclaunay distances 
(A^ e }=2.6 nm 



Longest Delaunay distances 



(A^ e )=25 



nm 



■ 1 1 1 1. ■ t ■ ■ 1 1 ■ , li i 1 1 1 > , u ■ , , , . . 



20 



40 



60 



Interprecipitate Distance (nm) 



FIG. 2: 3D IDD for the dataset in Fig. 1. (a) IDD of all 6,771 Delaunay lengths, with 
(A^? e )=16 nm. (b) Solid: IDD of nearest -neighbor distances, which is much sharper than when 
longer lengths are included «Af£ e )=2.6 nm). Hollow: IDD of the most-distant Delaunay neighbors, 
which is broader than and does not overlap with the shortest distances ((Ag^ e )=25 nm). 

a pair by a dislocation. The longest Delaunay distances provide an upper bound to the 
interactive distance. This is probably not physically important for plastic deformation, but 
may be relevant for other physical phenomena. The IDD for this case is broader, does not 
overlap the shortest distances, and has a mean value that is an order of magnitude larger 
«A^ e )=25 nm). 

In certain cases, it is not (Ag? e ) that is of interest, but rather the interplanar edge-to- 
edge distance, (Ag^ e ). This might, for instance, be a glide plane of a dislocation. This 2D 
distance can be calculated by imposing an additional constraint for COBYLA — that x and y 
values must fall on a particular plane. For comparison, (Ag? e ) and \\e- e ) can be calculated 
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FIG. 3: (a) Interplanar (2D) IDDs for the dataset in Fig. 1. Solid: IDD of slices, with 
<A^ e )=6.2 nm. Dashed: IDD from Fig. 2b scaled by Eq. 3 ((A2 D e }*=5.4 nm). (b) A^ e is 
weighted towards longer distances than \ 2 e D 



from one another by assuming precipitates are distributed on a cubic lattice [. 
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<«>=(viH <B>: (1) 

where is the volume fraction of precipitates. Assuming a square array of precipitates, 

<«} = yf-2)*; (2) 

where the mean planar radius, R, is equal to jU2(R), with U2 dependent on the PSD [2]. 
Values for u> 2 for the LSW distribution and for the case where all precipitates are the same 
size are given in Ref. [2]. Equating the 0s in Eqs. 1-2 leads to a cubic equation relating 
(Xl D e ) and <A^ e ). Solving for (\™ e ): 



tt -Auj 2 (Rf + JScul(R) (2(R) + (Xl D e )) 



<«> = m (3) 

In Fig. 3, the results of the two methods for extracting 2D nearest-neighbor IDDs from the 
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3D dataset are compared. \ 2 e _ e is calculated directly by imposing the additional constraint 
on COBYLA that the two points must lie in the same plane, which is radial to the analysis 
direction. The entire tip is sampled by taking 180 1° steps. This process samples the 
precipitates toward the center of the tip more than those toward the hull, but has fewer edge 
artifacts than parallel slices would have. \ 2r ? e * is calculated from what is displayed in Fig. 2b 
by applying Eq. 3 with (R) = 1.3 nm and o>2 = 1.046 [2] (the superscript * denotes this 
transformation). Despite the simple geometrical assumptions involved in deriving Eq. 3, the 
mean values are in reasonable agreement (with the "direct" method yielding (\ 2 ® e }=6.2 nm 
and with the conversion leading to (Ag? e )*=5.4 nm). Despite this similarity in the mean 
values, the distributions are different. The converted IDD is narrower and weighted towards 
shorter distances than the IDD that is directly calculated. 

This supports the conclusion reached in Ref . [1 1] , that A e _ e should be calculated directly 
with the same dimensionality as either the simulations they are compared with or the phys- 
ical models they are be used in. 2D techniques can only result in 2D IDDs, and will not 
yield accurate 3D IDDs. They may give a reasonable estimate of (X^ e ), although 3D ex- 
perimental data, as is gathered with the LEAP tomograph, allows both X^ e to be measured 
and Ag? e to be measured from planar slices taken from the 3D reconstruction. 

We are in the midst of applying this approach to calculate the strength of different alloys 
using analytical equations [2]. We are also using it to evaluate the statistical accuracy of 
simulated microstructures, which are used in a continuum dislocation dynamics simulation 
that calculates a stress-strain curve [31]. 
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